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Abstract —A key assumption made in this paper is that electric 
vehicle (EV) battery charging profiles are rectangular. This 
requires a specific and new formulation of the charging problem, 
involving discrete action sets for the EVs in particular. The 
considered cost function comprises of three components: the 
distribution transformer aging, the distribution energy losses, 
and a component inherent to the EV Itself (e.g., the battery 
charging monetary cost). Charging start times are determined 
by the proposed distributed algorithm, whose analysis is con¬ 
ducted by using game-theoretic tools such as ordinal potential 
games. Convergence of the proposed algorithm is shown to be 
guaranteed for some Important special cases. Remarkably, the 
performance loss w.r.t. the centralized solution is shown to be 
small. Simulations, based on realistic public data, allow one to 
gain further insights on the Issues of convergence and optimality 
loss and provide clear messages about the tradeoff associated 
with the presence of the three components in the considered cost 
function. While simulations show that the proposed charging 
policy performs quite similarly to existing (continuous) charging 
policies such as valley-filllng-type solutions when the non-EV 
demand forecast is perfect, they reveal an additional asset of 
rectangular profiles in presence of forecasting errors. 

Index Terms —EV charging - Energy scheduling - Transformer 
aging - Energy losses - Distributed algorithms - Game theory. 

I. Introduction 

The deployment of electric vehicles (EVs) at a large scale 
is envisioned to have a significant impact on the existing and 
future energy networks ||2l. In the present paper, the impact 
on the grid is assessed in terms of residential distribution 
network (DN) costs. To be more specific, the main goal 
pursued is to optimize EV charging schedules to minimize 
a cost resulting from a linear combination of the residential 
distribution transformer aging and the distribution energy 
losses. While the emphasis will be put on these two particular 
cost functions in this paper, all analytical results presented 
here apply to a large variety of costs; these results include the 
proposed problem formulation, the distributed charging algo¬ 
rithm, and its analysis. Among possible problems which can be 
considered we may mention voltage regulation IJI, harmonic 
distortion ||4|, and reactive power management Q. Concerning 
the component of the cost function which is referred as to 
the transformer component, the following has to be noted. 
Mathematically, it may correspond to any function of the past 
load levels; transformer aging is one possible instance and is 
the one made for the conducted numerical analysis. In practice, 
as explained in Ihl- lfT^ . optimizing a long-term criterion such 
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as the residential transformer lifetime becomes an important 
concern in the presence of EVs. Indeed, transformers might 
have to operate in a regime where aging is accelerated. In the 
European Union, about 5 millions of distribution transformers 
are used and about 70% of transformer failures are due to 
aging m. Concerning energy losses in the DN, they represent 
the most important fraction of power losses in the whole 
electricity network; according to Cl, in Erance two thirds 
of energy losses are due to the DN. Despite the importance 
of the problerrQ, the impact of EV charging on distribution 
transformer aging and energy losses has only been addressed 
in a relatively small number of papers. Among the relevant 
related works we may cite 1^ |[8l ifTOll ifTSll - lfTTll . The dominant 
approach adopted, which is well illustrated by 161, consists 
in exploiting a suitable model for the aging or energy losses, 
and assessing the impact of charging for simple scenarios', 
for instance, two possible scenarios are that all EVs start 
charging at a given time of the day (e.g., at 7 pm) or at random 
times. The algorithmic aspect is however not developed. This 
is precisely what the present work proposes. 

The algorithmic aspect of the charging problem has been 
tackled in the literature but, mostly, for minimizing the mon¬ 
etary cost the user has to pay for recharging his vehicle. In 
this respect llT8llfT9l constitute relevant works. Additionally, 
in the present paper the focus is on distributed charging in 
the sense that we assume the existence of several decision¬ 
makers (DMs) and each of them has only partial control of 
the variables which affect the cost or payoff function under 
consideration. Distributed charging is relevant in at least two 
key scenarios: the scenario in which charging policies have 
to be computed by a single entity (e.g., an aggregator ll20l or 
a transformer computing device) but for complexity issues it 
may be required to optimize the variables separately, forming 
a set of virtual DMs; the scenario where each EV controls its 
own charging policy, meaning that there are effectively several 
physical DMs instead of a single one. In the latter scenario, as 
far as advanced charging policies are concerned, DMs will be 
automata embarked on the EVs. Eor a naive policy such as just 
deciding when to plug the EV to the grid, the DMs might be 
the EV users themselves but this is not the standpoint adopted 
in this paper. 

It turns out that game theory is very well suitecJl to address 
the distributed charging problem whether players or DMs are 

* The problem directly concerns DN operators and car makers but may also 
concern the EV users since they might be charged in an indirect manner the 
extra costs induced by the impact of EV charging operations on the grid. 

^In particular, as explained further, the convergence and efficiency analyses 
for the proposed algorithm are conducted by introducing a charging game in 
which each player has his own payoff function. 
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persons or machines. In this respect, game-theoretic tools have 
been applied to smart grids quite recently (see e.g., ||2T]| for 
a survey); very useful contributions include 1^ . Il^ - ll25l . 
In ll22l . these tools are used for the DN frequency regulation 
problem in the context of the interaction between EVs and an 
aggregator. Therein, the authors show the usefulness of a well- 
chosen pricing policy to incite users to charge their vehicle in 
order to regulate the frequency of the electrical system. In Il20l . 
game theory is exploited to design a coordination mechanism 
for the wind power integration. References ll23l and ll24l use a 
similar methodology to study the more general problem of load 
balancing whereas ||25]| applies this methodology for micro¬ 
storage management in smart grids. 

Compared to the application-oriented works where game 
theory is used to optimize energy consumption at the user 
side (at home, by the EVs, etc.), the present work possesses 
several distinguishing features. Two of them are as follows. 
Eirst, to the best of our knowledge, our work is the first 
to propose a distributed charging algorithm which can not 
only minimize an individual cost, which is inherent to the 
EV and only depends on its actions (the battery charging 
monetary cost typically), but also the DN costs; the latter 
are given by the residential distribution transformer aging and 
energy losses over the DN. Second, we want to know to what 
extent using rectangular charging profiles is relevant both in 
terms of implementation and performance. The present work 
therefore adopts a complementary approach to existing works 
on charging algorithms which typically assume continuous 
charging power levels (see e.g., mum 122); the solutions 
used in the latter references are considered in the simulations 
section for comparison purposes. Although the focus of this 
paper is on the EV charging problem and DN costs, the 
obtained analytical and numerical results can be re-exploited 
for other problems in smart grids such as the problem of 
allocating or scheduling stored energy. 

The paper is structured as follows. Sec. |II] provides the 
model of the system and the network cost function components 
under consideration. Sec. HII] describes how the distributed 
EV charging problem is formulated. Sec. |IV] provides the 
description of the proposed distributed charging algorithm and 
its main properties. Sec. |V] allows one to assess numerically 
the performance of the developed algorithm. The paper is 
concluded by Sec. |Vl]. 

II. System modeling 

The goal of this section is to express the transformer aging 
and energy losses as a function of the sequence of the total 
demand power or load levels. To this end, the considered 
topology for the DN is as follows (see Eig. [U. It consists of 
one residential transformer to which two groups of devices 
are connected; a set of EVs and a set of other electrical 
devices. The latter are assumed to induce a power demand 
which is independent of the charging policies and therefore 
called exogenous demand. The corresponding load is denoted 
by Lf^°, which is a deterministic function of the time which 
is assumed to be slotted and whose index is denoted by f G 
T = {1, ...,T}; this function is always assumed to be known 


(except in the simulation part -Sec. |V]- where the influence of 
not forecasting it perfectly is assessed). The load induced by 
the EVs at time t will be denoted by Lf^ and the (average) 
total transformer load thus writes as Lt = Lf^° + Lf^. Note 
that our analysis concerns the impact of EV charging in terms 
of active power; refinements related to the reactive power are 
left as a possible extension. 



Eig. 1; Considered network topology. A typical exogenous 
(non-EV) demand profile Lf^°, is represented. 


A. Residential distribution transformer aging model 

The most influential parameter for the transformer aging is 
known to be the hot-spot (HS) temperature llZ7lll28]l . Indeed, 
the transformer isolation damage is directly related to the 
HS temperature: aging is accelerated (decelerated) when the 
HS temperature is above (resp. below) its nominal value. 
In general, the HS temperature depends on the history of 
transformer load levels and that of the ambient temperature 
levels. Even though the HS temperature at time t depends on 
the sequence of ambient temperature levels (see ED El for 
more detailsk we will not make this dependency explicit in 
our notation^ Indeed, this sequence cannot be controlled and 
in particular it does not depend on the power demand induced 
by the EVs. The HS temperature at time t is given by 

( 1 ) 

where i* = {Li,...,Lt) represents the sequence of total 
transformer load levels up to time f; the quantity Lt is 
related to the charging power levels in Sec. [Hll By default, 
no particular assumption will be made on the function F^^ 
throughout this paper. A few results will be provided for some 
special cases such as the memoryless case. Indeed, a special 
case of ([B is when the HS temperature does not depend on 
the whole history of load levels but only on the current load 
level; 

(^Lt). ( 2 ) 

Case o (resp. El) will be referred to as the case with (resp. 
without) thermal inertia. In E71 - E91 realistic models for the 
HS temperature evolution are provided, which allow one to 
have specific examples for F^^ and F^^. The correspond¬ 
ing specific expressions are only exploited in the numerical 
analysis. The transformer aging acceleration factor at time t, 
or aging for short, is assumed to be a function of the HS 

^The choice of a particular sequence of ambient temperature levels only 
intervenes in the numerical analysis. 
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temperature at time t and is denoted by At- Again, unless 
explicitly mentioned, no particular assumption will be made 
on the function At- A typical choice, which will be assumed 
in the numerical part, is as follows: 

= = (3) 

where a > 0, 6 < 0 are some constants (see CTlJll). 

Remark 1. All analytical results derived in this paper will 
only assume that At is a function of the past load levels 
(Li(s),Lt(s)). On the other hand, numerical results will 
we based on the particular choice coiTesponding to (3). 

B. Distribution network energy losses model 

Distribution energy losses mainly come from the trans¬ 
former and the lines between the transformer and the active 
electrical devices. At the transformer level, both voltage and 
frequency are assumed to be fixed and thus independent of 
the load level induced by the EVs. This is why no-load losses 
will not be considered here. As the purpose of this paper is 
not to obtain a model for energy losses which is as advanced 
as possible, we assume a simple scenario which is sufficient 
to show to what extent the mathematical problem formulation 
and charging scheme are impacted. Energy losses on time-slot 
t € T are assumed to express as 

J{Lt) = (i?ua„sfo + i?l.ne) {L^ + Lf^) ' (4) 

where i?tiansfo and i?iine are the transformer resistance and 
equivalent line resistance respectively; Lt = Lf^° + Lf^ is the 
total transformer load on time-slot t. As mentioned previously, 
the assumed model can be improved but the retained model 
has at least three attractive features: 1) provided that con¬ 
tinuous charging power levels are allowed, the minimization 
of Ji^t) w.r.t. corresponds to a valley¬ 

filling (VE) solution which is a well-known scheme ll26l . VE 
is performed over the sequence (L®™,..., L™°). Here this 
sequence corresponds to the exogenous demand prohle but 
it may also represent a sequence of prices (see e.g., mi); 2) 
the assumed energy losses model allows us to leave the load 
flow problem as a separate problem which might be handled 
with through an extension of our work lf30l : 3) the model is 
relevant when line energy losses are dominated by transformer 
energy losses. 

III. Formulation of the EV charging problem 

This section aims at formulating in a distributed manner (as 
motivated in Sec. |I] the EV charging problem. The problem is 
said to be distributed because the variables which affect the 
payoff (or cost) functions of interest are not controlled jointly 
but separately. Since the EV charging profiles are imposed to 
be rectangular, these variables correspond to the charging start 
times of the different EVs. The number of consecutive time 
instances or time-slots required to have the battery charged 
or to reach a required state of charge (SoC) for the next trip 
of EV i is denoted by Ci while the effective charging start 
time for EV i is denoted by Si. The individual payoff function 


which has to be maximized for EV i G I, X = is 

assumed to have the following form: 

= -/, {gf^{s)+gf^{si)) ( 5 ) 

where s = (si,..., sj) and ft is the individual pricing function 
for user i, which is assumed to be strictly increasing; for 
instance, it may translate the technological costs induced by 
charging into a monetary cost. The function (s) represents 
the cost associated with the DN. This cost is chosen to be 
a linear combination of the transformer aging and energy 
losses; its exact expression is provided a little further. The 
function gf^ can be any single-variable function of Si. It 
is an individual cost which only concerns EV i. It may 
model the impact of the start time on EV i battery aging, the 
individual electricity fare for user i, or its preference in terms 
of availability. For example, if the sequence of prices for EV 
i over the time period of interest is denoted by (71^1 ,..., t) 
(i.e., the price is a function of the time only), then a suitable 
choice for gf^ might be 

Si+Ci-l 

gf'^isi) =13 7r,,t (6) 

t—Si 

where /3 > 0 is a weight which allows the tradeoff between 
individual EV preferences and DN costs to be tuned. The 
price model can even be more complicated mathematically 
than what is assumed in © e.g., by assuming that t is a 
function of the current number of EVs charging at time t i.e., 
T^i,t = 4>i,titit) (see IIT9I for more details on this model). One 
of the main properties (namely, potentiality) of the charging 
game studied in Sec. |IV]is retained. 

To explicit gP^ as a function of s, some notations have to be 
introduced. Let n* and rit respectively denote the numbers of 
EVs starting to charge and charging at time t. These quantities 
are related to s = (si,..., s/) by 

/ I Ci 

^t{s) = Y 1^=*] = X! X! (7) 

2—1 2—1 t' — l 

where 1[ ] is the indicator function. To evaluate the impact 
of the charging policies on the DN, the key quantity is the total 
transformer load or consumed power. The T—dimensional 
sequence of total load levels is now also denoted by L'^{s) = 
{Li{s),L t{s)) and the total load on time-slot t is given 
by 

Lt{s) = LT° + Tp^(s) = LT° + Pnt{s) (8) 

where P is the common charging power of all the EVs. Using 
the introduced notations, it is now possible to express gf^'- 

g?^{^)= Y oiA{L*{s)) + {l-a)J{Lt{s)) (9) 

teWi(si) 

where 0 < a < 1 is the weight given to transformer aging 
relatively to energy losses and Wi{si) is a discrete set which 
represents the time window over which the EV i is considered 
to be influential on the DN cost. We will dedicate more 
attention to two special cases of practical interest for Wi(si). 
The first case is when yVi(si) = {si, • • ■ , -f — 1} = Wi 
which means that each EV individual payoff is only related 


TRANSACTIONS ON SMART GRID 


4 


to the period of time over which the corresponding EV is 
active. The second case corresponds to yVi(si) = T. The 
parameter a can be seen as a simple way of tuning the tradeoff 
between a short-term cost (energy losses) and a long-term cost 
(transformer aging). 

Now let us specify the action set for each EV. The arrival 
and departure time of EV i € I are denoted by Oi G T and 
di &T respectively. As it is assumed that an EV has to charge 
its battery within the total time window, the action set for EV 
i is chosen to be 

Si = {ai, Qi 1,di — Ci + 1} . (10) 

The EV action profile s therefore lies in 5 = 0^=1 Si', the 
standard notation = (si,..., Si_i, Si+i,.., s/), / > 2, will 
be used for referring to the reduced action profile in which 
user i’s action is removed. A special case of interest is when 
all users have the same charging constraint, that is Vi, at = 
a, di = d,Ci = C. This case will be said to be symmetric. 
In the symmetric case, it will be assumed, without loss of 
generality, that a = 1 and d = T. 

Remark 2. Maximizing the sum of weighted individual costs 
provides a Pareto optimal point when the cost region is convex 
ED. By changing the weights, one moves along the Pareto 
frontier of the cost region, which represents the best which 
can be done in terms of tradeoff between aging and energy 
losses. This therefore gives a motivation for considering the 
linear combination. A natural question is about whether the 
cost region is convex. In fact, the feasible average cost region 
is convex whenever the cost is averaged over a large number 
of stages (which is the number of days here). This follows by 
a time-sharing argument; if there exist two charging strategies 
which achieve each a given pair of individual costs, then any 
convex combination of these pairs can be achieved by using 
the two strategies with the appropriate fraction of the time. 

Remark 3. As already mentioned, charging profiles are 
assumed to be rectangular, which is why the charging power 
can only have two possible levels namely either 0 or P. 
Among the motivations for considering rectangular charging 
profiles we may mention the following; 1) An important 
argument is that rectangular profiles are currently being used 
for existing EVs (e.g., the EVs built by the Erench car 
maker RENAULT) and not only in papers; 2) for a given 
charging start time, charging at full power without interruption 
minimizes the delay to charge; 3) one technological reason is 
that they allow one to manage the EV battery aging. Battery 
aging seems to be accelerated when the charging operation 
comprises interruptions EH; 4) as shown further in Sec. |V] 
rectangular profiles are also fully relevant for cost functions 
with memory such as the transformer aging. The charging start 
time turns out to be a very influential variable, confirming 
the observations made in Col; 5) more specifically, from an 
optimal control theory perspective, rectangular profiles may be 
optimal. This happens for instance when the state (i.e., the hot¬ 
spot temperature) is monotonically increasing with the control 
(i.e., the charging power). If, for the time window of interest, 
the transformer temperature can only increase, it is optimal 
to delay the transformer heating. The optimal solution is then 


to start charging as late as possible i.e., to charge at maximal 
power at the end of the considered time window and charge at 
zero power before. The duration of the corresponding charging 
profile, which is rectangular, is given by the desired final state 
of charge; 6) the charging start time will be seen to be less 
sensitive to forecasting errors on the exogenous demand than 
VE-type solutions; 7) profiles without interruption are even 
required in some important scenarios encountered with home 
energy management E21E31- 

IV. A NEW DISTRIBUTED CHARGING ALGORITHM 

A. Motivations 

As explained in Sec. HI considering distributed charging 
policies is relevant in at least two key scenarios. Assume a 
scenario (say Scenario c) in which the charging policies are 
computed by a single decision-making entity (e.g., an aggrega¬ 
tor or a transformer computing device) and the maximization 
of a quantity such as the sum-payoff function X]i=i 
pursued. Einding an optimal solution may largely exceed the 
available computational capacity; an exhaustive search would 
roughly involve tests. Therefore, even if there is one single 
decision-making entity, it may be required to optimize the 
variables of s separately. Now, if we assume a scenario (say 
Scenario d) in which each EV controls its own charging policy, 
there are I physical DMs and the problem is distributed by 
nature. In both scenarios, the total computational complexity 
of the proposed algorithm will be seen to be typically linear in 
the product T x I, showing the dramatic reduction in terms of 
complexity allowed by the used distributed implementation. 

The proposed distributed algorithm is based on a procedure 
which is called the sequential best-response dynamics (BRD) 
in game theory literature (see e.g., ESESl). One of the strong 
motivations for selecting such a procedure for the problem un¬ 
der investigation is that convergence of the associated iterative 
algorithm can be guaranteed with overwhelming probability. 
Additionally, convergence is very fast, which is useful both in 
Scenario c to avoid unnecessary computations and in Scenario 
d to manage the amount of signaling between the EVs and the 
aggregator. Other arguments in favor of using the BRD will 
be provided further. 

To clearly indicate that no strategic assumption such as 
rationality or complete information is required on the DMs 
which implement the BRD, the description of the algorithm 
has been separated from its analysis. The analysis relies on 
the use of game-theoretic tools such as the powerful notion of 
potentiality, which guarantees the existence of a Nash point in 
a game and the convergence of the BRD to a Nash point. 
It turns out that the charging game under consideration is 
effectively potential and therefore makes the BRD a good 
candidate for computing the charging start instants. 

B. Description of the algorithm 

The proposed algorithm to determine the vector of start 
times (si,...,s/) is an iterative algorithm which is inspired 
from the sequential BRD. The algorithm is performed offline, 
which means that the decisions which intervene in the algo¬ 
rithm are intentions but not decisions which have effectively 
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been taken; only the decisions obtained after convergence 
will be effective and implemented online. Once the charging 
instants are computed, the EV can effectively charge their 
battery according to the schedule determined. In its most used 
form, the BRD operates sequentially such that DMs update 
their strategies in a round-robin manner. Within round m (with 
m > 1) the action chosen by DM i (which can be virtual 
or physical depending on the assumed scenario) is computed 
as (dUi- The proposed procedure is translated in pseudo-code 
through Algorithm 1. 


Algorithm 1: The proposed distributed EV charging algo- 
rithm._ 

Initialize the round index as m = 0. Initialize the vector 

of charging start times as 

while II— s(’”“^^|| > 5 and m < M do 

Outer loop. Iterate on the round robin phase index: 
m = m + 1. Set i = 0. 

Inner loop. Iterate on the DM index: i = i + 1. Do: 

(m) _ / (m) (m) 

s- ^ € arg maxiii(5i % S 2 

Si^Si 

( 11 ) 

where stands for action of DM i in the round 
robin phase m. Stop when i = I and go to Outer 
loop, 
end 


Comments on Algorithm 1. 

• In ([IT]l, when the argmax set is not a singleton, is 
randomly drawn among the maximum points. 

• The quantity <5 > 0 in Algorithm 1 corresponds to the 
accuracy level wanted for the stopping criteria in terms of 
convergence. 

• To update the charging power levels m times, mxl iterations 
are required. 

• The order in which DMs update their action does not matter 
to obtain convergence (see e.g., ETI). However, simulations 
which are not provided here indicate that some gain in terms 
of convergence time can be obtained by choosing the order 
properly. A good rule seems to be to start updating at iteration 
m the EV decisions in an increasing order in terms of start 
times as obtained per iteration m — 1, which makes the order 
iteration-dependent. 

• The knowledge required to implement Algorithm 1 is 

scenario-dependent. In Scenario c in which each decision is 
computed by a single entity (the transformer typically), the 
vector of effective charging instants can be computed from 
its initial value the (forecasted) sequence of exogenous 
loads ...,L^°), and the parameters which intervene in 

the payoff functions; the latter include in particular the EV 
mobility data {ai)i^x, {di)i^x, and (C'i)igx- In Scenario d 
in which the EVs themselves update their decision, messages 
have necessarily to be exchanged between the transformer and 
the EVs. A possible communication protocol is as follows. 


Without knowing anything about the exogenous demand 
and the one associated with the other EVs, EV automaton i 
chooses a start time say and reports this to the transformer. 
The latter aggregates the received signals and replies to the 
EVs by sending them the predicted sequence of total load 
levels = Lf^° -f Put^), t € T. Therefore, EV 1 
updates its intended start time as € argmax^j rti(si, 
and reports this change to the transformer: the latter updates 
the aggregate signal into = Lf^° + Pnt{s^^^), t € T, 
where ..., This signal is sent to all 

the EVs but only EV 2 is able to update its intended start 
time. When all EVs have updated their start time at least 
once, a new updating round can start. In practice, it might 
happen that much less knowledge is available. Eor instance, if 
only the knowledge of the forecast exogenous demand 
is available then it is always possible to apply Algorithm 1 
on the cost function 0 with a = 0. If the corresponding 
charging scheme is used, it will induce, in general, a certain 
loss of optimality. 

• A variation of Algorithm 1 can be obtained by updating the 
charging policies simultaneously. The main reason why we 
have not considered the parallel version is that it is known 
that there is no general analytical result for guaranteeing con¬ 
vergence ll^ . When converging, the parallel implementation 
is faster but since start times are computed offline, convergence 
time may be seen as a secondary feature. 

C. Convergence analysis 

One of the powerful links between distributed optimization 
and game theory is that scenarios involving several individual 
optimizers or DMs which update their decisions over time may 
converge to a Nash equilibrium (NE) of a certain game. This is 
one of the reasons why we now define the game of interest, that 
we will refer to as the charging game. The main purpose of this 
section is to show that, under additional realistic assumptions, 
this game is an ordinal potential game (OPG) IITSl . The 
intuition behind this is that the considered cost functions all 
depend on the EV charging loads through the sum-load. This 
type of structures, which is present in games which are called 
aggregate games (a congestion game ll39l is a special instance 
of them), may lead to the existence of a potential function. 
Although exact potentiality is typically not available here, it 
turns out that ordinal potentiality is available under some mild 
conditions. In particular, the latter property guarantees the 
convergence of the distributed charging algorithm proposed in 
Sec. |IV] A game under strategic-form is given by an ordered 
triplet which respectively comprises the set of DMs or players, 
their strategy sets, and their payoff functions (see e.g., ED). 
The charging game of interest is defined as follows. 

Definition 1 (Charging game). The charging game is the 
triplet Q — (l, ^dose elements are defined 

by 0 and ( 17^ . 

Interestingly, this game is an OPG for typical scenarios 
encountered in practice. An OPG is defined as follows lf38l . 

Definition 2. [Ordinal Potential Game] A game whose payoff 
functions are is an OPG if there exists a function $ 
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such that \/i G I, Vs = (sj, S-i), Vs' G 

Mi(s', S_i) > Mi(si, S_i) <^=> $(s', S-i) > <i)(si, S_i) . 

Since the function $ does not depend on the player index, 
the game analysis amounts, to a large extent, to analyzing an 
optimization problem. This attractive property is available at 
least in the scenarios defined by the next proposition. 

Proposition IV.l (Potential property of the charging game). 
If one of the following conditions is met, then Q is an 
OPG: (a) Vi G I, Wi(si) = W, where W is any discrete 
set which is independent of the player index i and Si; (b) 
Vi G I, yVi{si) = {si,--- ,Si + Ci — 1} and there is 
no thermal inertia i.e., 0^® is given by m- Two potential 
function candidates are respectively given by: (a) $( 3 ) (s) = 

- (^)) + (1 - (s)) - 

tew iei 

ntis) 

^(b) (*) = - XI i^tivt)) + (1 - a)J {Lt(vt)) - 

teT vt=o 

X Ltivt) = Lf^° + Pvt. 

iei 

The proof is provided in App. The two scenarios in which 
the game is potential are clearly of practical interest. Assuming 
all the EVs to have a common time window Wi(si) = W for 
charging (Scenario (a)) means that the total cost associated 
with the DN is accounted for by all the EVs. This is clearly 
the most interesting scenario for the aggregator. However, from 
the user’s standpoint, this may be considered as unfair or 
not acceptable. In such a case, it is more realistic that user 
i be only charged a cost which corresponds to the period over 
which his battery is effectively recharged namely, in the time 
interval defined by Wi{si) = {si, Si + 1,..., Si-\-Ci — 1} (note 
that an alternative way of individualizing the cost would be 
to use a cost sharing policy such as in which would lead 
us to use a common window W but weighting the total cost 
by for user i). This general scenario is mathematically 

more involving than Scenario (a). It turns out that it becomes 
quite simple to be analyzed for transformers with low thermal 
inertia, which leads to Scenario (b). Obviously, when energy 
losses represent the dominant cost (a ^ 0), the game is always 
potential. All these comments lead us to the next proposition 
which is the key result of this section. 

Proposition IV.l (Sufficient conditions for convergence). Al- 
gorithm\J]converges if one of the following conditions is met: 
(a); (bj; a = 0. 

As proved in 1^ . the sequential BRD converges in ordi¬ 
nal potential games. The proof therefore follows. When the 
charging game does not meet none of the conditions above, a 
deeper analysis has to be conducted. In this paper, our choice 
is to exploit Monte-Carlo simulations to provide additional 
insights on the convergence issue. In Sec. [V) the empirical 
convergence probability is assessed for other typical scenarios 
and is shown to be high. When converging, the proposed 
algorithm converges to an NE of the charging game Q. The 
motivation of the next section is to analyze the existence, 
uniqueness, and global efficiency of the convergence point(s) 


which are NE. 

D. Equilibrium analysis 

A pure NE of is a point which meets a certain condition 
of stability l40l . Eormally, it is defined as follows. 

Definition 3 (Pure NE). The action profile s* = (s*,..., s}) G 
S is a pure NE if \/i,\/si, Ui{si,s*_^) < Ui{s*). 

The existence of a pure NE in Q is ensured under the 
conditions assumed in the following proposition. 

Proposition IV.3 (Existence of a pure NE). In Scenarios (a), 
(b), or a = 0, the game Q has at least one pure NE. 

The proof of this result follows from the fact that in any 
of the scenarios of Prop. II V. 1 1 Q is an OPG and after ll^ 
the existence of a pure NE is guaranteed. On the other hand, 
and this is common in discrete games, uniqueness does not 
hold. To prove this set T = 5, / = 3, = C = 2, = 
a = 1, d, = d = 5, = (1, 2,3, 2,1), and P = 1. Eor a 

small transformer inertia i.e., a small (see Sec. |V]l, it can 
be checked that s* = (1,1,4), s** = (1,4,1), and s*** = 
(4,1,1) are NE, showing that uniqueness is not guaranteed 
in general. Since NE uniqueness is not guaranteed in general, 
measuring the efficiency of the worst NE is important. The 
most usual way of assessing the impact of decentralization 
in a non-cooperative game has been formalized in m by 
defining the notion of price of anarchy. Rather, we will slightly 
modify the latter notion as the price of decentralization (PoD), 
which we define below; the merit of the proposed definition 
is just that the price is effectively zero when a distributed 
algorithm or procedure leads to an equilibrium point which 
performs as well as the centralized solution in terms of sum- 
payoff w = 

Definition 4 (PoD). The PoD of Q is defined by 

max w{s) 

PoD = 1 - -— (12) 

min w(s) 

seSNE 

where 5^® is the set of NE of the game. 

It can be seen that 0 < PoD < 1 and the larger the 
PoD, the larger the loss due to decentralization. It is generally 
difficult to express the above quantity as a function of the game 
parameters ll^ . This explains why these quantities are often 
more relevant from the numerical point of view. Nonetheless, 
it is possible to characterize it in some special cases. One of 
the cases where PoD can be characterized is the limit case 
of a large number of EVs, that is / ^ 00 , having the same 
charging constraint Oi = a, di = d and Ci = C, and the 
transformer has no thermal inertia. In this asymptotic regime, 
^ Xt G K represents the proportion of EVs charging at 
time t and the analysis of the game Q amounts to analyzing the 
so called non-atomic counterpart of Q. In the latter game, 
the set of players is continuous and given by = [0,1]. 
The action set of the EVs, S, is defined as in (fTOl i. In the 
regime of large numbers of EVs, the transformer load becomes 
L{x) = -\-px; the parameter p is introduced in order for 
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the exogenous demand to scale with I. Indeed, after ([8]), when 
/ —+ 00 , if kept fixed, the exogenous demand tends to 
vanish in comparison to the load induced by the EVs. This is 
the reason why we introduce the parameter p (instead of P). 
The obtained non-atomic charging game can be proved to be 
an OPG and the following result concerning efficiency can be 
obtained. 

Proposition rV.4 (PoD in the non-atomic case (I oo)). 
Assume that: ai = a, di = d, Ci = C, and Wi(si) = T; 
gf^ = 0; is a non-increasing (resp. non-decreasing) 

function of t on {1, • • ■ , C} ( resp. {T — C + I,-- - ,T}); 
Scenario (b) or a = 0 is considered. Then we have that 
PoD = 0. 

The proof of this result is provided in App. |B] This result 
has the merit to exhibit a scenario where decentralizing the 
charging decisions induces no cost in terms of global optimal¬ 
ity for the sum-payoff. Note that, in particular, if the exogenous 
demand is either constant or negligible w.r.t. the demand 
associated with the set of EVs, the above assumption holds and 
there is therefore no efficiency loss due to decentralization. 

Remark 2. The fact that there exist multiple Nash equilibria 
might be seen as a crucial point since one does not not know to 
which point Algorithm 1 will converge. However, Prop. IV.4 
and the numerical results provided in the next section show 
that the PoD is typically small. This means that the worst Nash 
equilibrium and the best Nash equilibrium necessarily perform 
similarly, showing that equilibrium selection is not a crucial 
issue for the problem under consideration. 


A6iPl = 55°C; = 23°C; q = V, r = f 6^^ = 98°C. 

Energy losses are evaluated by choosing i?transfo = ^line = 
0.03 n. The load induced by one EV is P = 3 kW. Eor each 
day, charging operations have to take place within the time 
window from 5 pm (day number j) to 8 am (day number 
j -f 1), which corresponds to T = {1, 2,..., 30}. Concerning 
the EV mobility data, two scenarios will be considered. In 
Scenario (s), all EVs need = C = 16 time-slots of 30 
minutes each to completely recharge their 24—kWh battery 
and Oi = a = 1 (5 pm), di = d = 30 (8 am); this scenario 
can be seen as the worst case. In Scenario (t), the mobility 
data are deduced from statistics taken from the Erench survey 
ENTD 2008 available in ll44l : Oi, di, and Ci are taken to be the 
closest integers of realizations of Gaussian random variables 
a* ~ A/'(4,1.5), d^ ~ A/'(29,0.75) and C* ~ A/'(5.99,1.14). 
By default. Scenario (s) will be assumed. Einally, unless 
specified otherwise, the cost functions in (|5]l are defined by: 
for all i, gf'^ = 0, which allows us to isolate the effects of the 
exogenous demand; Wgsi) = Si-f 1,..., -I-— 1}, 

and fi = Id. The plug-and-charge (PaC) policy is obtained 
by assuming that EVs start charging as soon as they plug to 
the grid according to the data of ll44l . To assess the impact 
of not being able to forecast the non-EV demand perfectly 
we assume that the available non-EV demand profile is given 
by Lf^° = Lf^° Z, where Z ~ j\/"(0, We define the 

forecasting signal-to-noise ratio by 

/ , , Tday \ 

FSNR = lOlogio — X (dB) (14) 

y^^day t^i I 


V. Numerical analysis 


We first provide the general simulation setup assumed by 
default while particular choices will be specified in the figure 
captions. Data corresponding to non-EV demand (or exoge¬ 
nous) profiles and the ambient temperature are taken from the 
ERDE Erench DN operator data basis. They concern Erance for 
the year 2012 and can be found in ll42l . Il43l . Unless specified 
otherwise, simulations are performed over a year; the chosen 
time unit corresponds to 30 min. We consider a 20 kV/410 V 
transformer whose apparent power is 100 kVA and nominal 
(active) power is 90 kW. The transformer HS temperature 
evolution law is assumed to follow the ANSI/IEEE linearized 
Clause 7 top-oil-rise model, which is described in ll29l . The 
transformer lifetime is inversely proportional to the average 
aging: 

/Uear \ 

lifetime = 40 x Tyear x j ^ At | (years) (13) 


where the non-EV or exogenous demand is normalized such 
that without EV lifetime = 40 years; here Tyear = 366 x T. 
The instantaneous aging At (see Q) is computed by choosing 
a = 0.12°C“^, b = —11. The function is not described 
here but can be found in ll29l . To make the simulations 
reproducible we provide the values of the different parameters 
of At = 0.5 h; T° = 2.5h (thermal inertia) for all simu¬ 
lations concerning the transformer, unless specified otherwise 
(in which case we have that T° = 0.5h) ; 7 = 0.83; R = 5.5; 


where Tday = 48. 

Numerical convergence analysis. In Sec. IIV-CI we have 
provided sufficient conditions under which convergence is 
guaranteed. Here, we consider a scenario in which these 
conditions are not met (in fact the worst, with a = 1 ) and 
assess the probability of convergence of Algorithm 1. Eor 
Scenario (s), Eig. |2] represents the empirical probability of 
convergence against the number of EVs for the 366 exogenous 
demand profiles from Il42l and for 10 000 draws from a 
Gaussian random vector; the covariance matrix of the latter 
is taken to be x I{txT}^ where a = 26kW is estimated by 
using the data from ll42l . All simulations performed showed 
that only a few iterations are needed to obtain convergence, 
which is a quite typical behavior for sequential BRD-type 
iterative procedures 13^ . 

To elaborate further on this point, we provide Eig. |3] It rep¬ 
resents the total number of iterations needed for convergence 
of Algorithm 1 as a function of the number of EVs when 
Wi(si) = T. The middle curve corresponds to an average over 
the 366 days of 2012 while the two others are associated with 
a chosen confidence interval of 68 %. The total convergence 
time is observed to scale very well with the number of EVs 
indicating no scalability issues for convergence. Note that 
the scaling law seems almost piecewise linear here but many 
other simulations performed for more diverse scenarios show 
it might be "less linear" but always involves smooth variations. 

To conclude on the convergence issue, we consider the vari¬ 
ations of the individual payoff functions of the EVs and that 
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Fig. 2: Even if the assumption of Prop. IIV.2I is not met. 
Algorithm 1 still converges with high probability. When the 
number of EVs (/) is low or large, the (empirical) probability 
of convergence tends to 1. When the number of EVs takes 
intermediate values, the probability is typically above 90%. 



Eig. 3; Total number of iterations of Algorithm 1 against the 
number of EVs (I). The dashed curves represent the 68% 
confidence intervals. It is seen that Algorithm 1 is scalable 
in terms of convergence time. Here convergence time scales 
well with I since variations are smooth. 

of a potential function over iterations; this is the purpose of 
Eig. 13 Eor / = 10, Vj, = 1, di = ... = dg = 30, dig = 24, 
the night of January 1, 2012 ll42l . and now a = 0 (only 
energy losses are taken into account), the figure shows that the 
potential always decreases over iterations, which illustrates the 
fact that it is a Lyapunov function of the considered dynamical 
procedure. As a very positive result, it is seen here that only 
one update per EV is needed to reach convergence. This very 
fast convergence behavior is typical when a memoryless cost 
is considered which is the case here with a = 0. 

HS temperature evolution. Eor Scenario (s), Eig.|5]represents 
the HS temperature against time for the worst day of the year 
2012 in Erance in terms of HS temperature peak (top curve) 



Eig. 4: Eor 10 EVs, only 10 iterations are needed to converge, 
which means that only one update per EV is sufficient to reach 
convergence, showing a very fast convergence. The figure also 
illustrates the fact that the potential is a Lyapunov function of 
the considered dynamical procedure. 

and the best day (bottom curve) in three scenarios: without EV, 
with Algorithm 1, and with the PaC policy. The selected worst 
day (which is a winter day) shows that the HS temperature 
takes excessive values when the PaC policy is used. The 
temperature peak is increased by about 80°C w.r.t. to the case 
without EV (see na for physical justifications on the value 
of this excess). On the other hand, implementing Algorithm 
1 does not induce any increase for the peak and roughly 
tends to minimize HS temperature variations. Since the aging 
acceleration factor At is exponential in the HS temperature, 
the PaC policy has a dramatic effect in terms of transformer 
lifetime. Note that this decrease might be made much slower 
by choosing a transformer with a larger admissible active 
power (e.g., 120 kW instead of 90 kW) but here, our goal 
is to see the impact of EVs on residential transformers which 
are already deployed and have therefore been sized to operate 
without EV. 

Performance comparison analysis and influence of the fore¬ 
casting noise. Here, we set a to one. The transformer is 
assumed to be chosen to be able to operate for 40 years without 
EV. Eor Scenario (t), Eig.|6]represents the transformer lifetime 
against the number of EVs when the exogenous demand 
forecast is perfect (see the three curves in dotted lines) and 
when FSNR = 4 dB (see the three curves in solid lines). The 
top curve is an horizontal line which corresponds to the case 
without EV; the bottom curve corresponds to the PaC policy. 
The three non-trivial charging policies under consideration are 
the one corresponding to Algorithm 1, that of Gan et al ifT^ 
and Shinwari et al 12^ . The Gan et al policy corresponds to 
the convergence point of an iterative algorithm which aims at 
minimizing a cost which results of two terms: if maximized 
alone and assuming convergence, the first term would lead to 
a VP solution; a second term whose role is to stabilize the 
parallel implementation-based iterative algorithm. The weight 
put on the latter penalty term (0.5) is tuned optimally for Pig. 
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Fig. 5: Hot-spot temperature against time without EV (for 
two extreme days in 2012), with Algorithm 1, and with 
the PaC policy. Scheduling charging start times according 
to Algorithm 1 instead of plugging-and-charging allows the 
temperature variations to be minimized and does not induce 
any excess for the peak. It is also seen that scheduling charging 
needs properly is much more important during winter time (in 
France). 

|6l if this weight is not tuned properly, the implementation of 
Gan et al may lead to significant performance losses ll46l . For 
the Shinwari et al policy, the energy need of EV i is spread 
by filling the "holes" of the exogenous demand. For each EV, 
a proportion of the energy needed is allocated to a given time- 
slot proportionally to ^ with St = + max* 1/?’^°, 

and the remainder is uniformly allocated. Fig. shows that 
the PaC policy is seen to be non-acceptable, showing the 
imperious need for advanced charging schemes. It is seen 
that Algorithm 1, which is based on rectangular charging 
profiles, performs as well as the continuous power level-based 
scheme of Gan et al. Both schemes yield a relatively small 
decrease of transformer lifetime in Scenario (t). This holds 
under the assumption of perfect forecasting for the non-EV 
demand profile. However, under the more realistic assumption 
of imperfect forecasting (FSNR = 4 dB), as seen in Fig. |6] 
transformer lifetime can severely be degraded as the number 
of EVs increases for the two latter schemes. Algorithm 1, 
which is based on rectangular charging prohles, is seen to be 
much more robust against noise on the exogenous demand 
than VF solutions. This observation is clearly confirmed by 
Fig. |7] which is also obtained in Scenario (t) and represents the 
transformer lifetime against FSNR for three charging schemes. 
Rectangular profiles have the advantage to be less sensitive to 
amplitude errors than the VF solution since the sole parameter 
to be tuned is the charging start time. We have checked that this 
message is not mobility data-dependent. Indeed, simulations 
which are not provided here and exploit the (US) mobility data 
of in and Bsl confirm that rectangular charging profiles are 
robust against forecasting errors. This provides a very strong 



Number of electric vehicles 


Fig. 6; Transformer lifetime against the number of EVs (/). 
The plug-and-charge policy is seen to be non-acceptable. 
Algorithm 1 is seen to perform as well as existing valley-filling 
type solutions under perfect forecasting of the non-EV demand 
and outperforms these solutions under imperfect forecasting. 



Fig. 7: Transformer lifetime (years) against the forecasting 
signal-to-noise ratio (dB) for / = 10 electric vehicles; a = 1 
and mobility data of Scenario (t), namely those from Il44l . are 
assumed. 

argument in favor of using rectangular profiles. A parallel with 
the problem of noise robustness for high-order modulations in 
digital communications can be drawn and an optimal power 
level might be identified, which is the purpose of the next 
paragraph. 

Existence of an optimal charging power level for rectangu¬ 
lar profiles. While it is clear that the larger the charging power 
the lower the time to charge, having a high charging power 
can be suboptimal in the presence of forecasting noise. This 
is what Fig. [8] shows in Scenario (s). It depicts the optimal 
charging power in terms of transformer lifetime at the NE 
obtained with Algorithm 1 as a function of the number of 
EVs, the minimum (resp. maximum) charging power being set 
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Fig. 8: In the presence of forecasting noise, charging at the 
highest charging power possible is suboptimal in terms of 
transformer lifetime. Rather, a non-trivial optimal charging 
power level can be determined. 

to 2.2 kW (resp. 24 kW). At one extreme, in the absence of 
forecasting noise, the optimal power level is always 24 kW, 
which corresponds to two 30—min time-slots. At the other 
extreme, when the noise level is large, the best power level 
would correspond to charge over the entire period namely at 
= 1.6 kW if it was allowed and we find the minimal 
value of 2.2 kW. For typical forecasting noise levels Il49l . 
simulations reveal non-trivial optimal charging power levels, 
as illustrated by Fig. |8] 

Tradeoff between the transformer lifetime and energy losses. 
Fig. |9] represents the Pareto frontier of the feasible cost region 
for the first day of 2012 when the proposed algorithm is used; 
the x-axis corresponds to the normalized transformer aging 
while the y-axis corresponds to normalized energy losses. The 
curves are obtained by varying a from 0 to 1 and considering 
two different scenarios in terms of thermal inertia: = 0.5 h 

(top curve) and T° = 2.5 h (bottom curve). The conclusion is 
that it is preferable to design charging policies which minimize 
the transformer aging, i.e., to set a = 1. The loss of optimality 
in terms of energy losses will be rather small by using these 
policies. It has been observed that changing the charging start 
time by one hour or two does only affect energy losses in a 
marginal way; Tab.|Il which provides normalized energy losses 
for different charging schemes and number of EVs, confirms 
this. However, because of thermal inertia, changing the start 
time by one hour or more has a significant impact on aging. 
This is due to the fact that typical exogenous demand profiles 
comprise a valley in the night, explaining the results of Fig. 

in 

Tradeoff between transformer aging/energy losses and 
charging monetary cost. The purpose of Fig. [TO] and Fig. fTTIis 
to assess what is lost in terms of charging monetary cost when 
pursuing transformer aging or energy losses minimization. 
Indeed, these figures depict the charging monetary cost against 
the number of EVs when: an EV aims at minimizing the 



Normalized transformer aging 

Eig. 9: Charging start times have an important impact in terms 
of transformer aging while they have much less influence on 
energy losses. 


J 

5 

10 

20 

30 

Plug-and-Charge 

1.14 

1.30 

1.70 

2.18 

Algorithm 1 

1.09 

1.21 

1.50 

1.86 

Gan et. al. 0181 

1.09 

1.20 

1.49 

1.84 

Shinwari et. al. 0261 

1.10 

1.22 

1.50 

1.85 


TABLE I: Normalized energy losses J (they are normalized 
relatively to the case without EV) against the number of 
electric vehicles (/) for four charging schemes. 

charging monetary cost (we force to be zero in (|5]l and 
only exploit the function gf^f, or it aims at minimizing the 
transformer aging (i.e., when a = 1 in (|9l) and gf^ = const.); 
or it aims at minimizing energy losses (i.e., when a = 0 in (|9]l 
and gf^ = const.). This is done for two choices of electricity 
fares: the Erench on/off peak fare; the market price (Epex Spot 
prices in Erance ll50l ). It can be seen that choosing a good 
charging scheme in terms of aging or energy losses leads to 
a charging monetary cost which is reasonably close to the 
one which is obtained by minimizing the monetary cost. The 
explanation for this is as follows. If the local demand (e.g., at 
the residential scale) is correlated to the global demand (at the 
scale of the country), a good tradeoff between technological 
costs and the monetary cost can be found. Indeed, when 
looking at the non-EV demand and the electricity price profile, 
it is seen in particular that there is an "evening peak" both 
for the demand and the price. When using other data like 
the ERGOT (Electricity Reliability Council of Texas, ifSTI ) 
data, similar conclusions can be drawn. Now, if the non-EV 
demand and prices are not correlated at all, the considered 
cost function has to combine three terms and the value of 
the weight /3 (relatively to a) will be very influential on the 
charging profiles obtained with Algorithm 1. Corresponding 
simulations might then be provided but as our initial goal was 
to put the emphasis on the distribution network costs and not 
on the monetary aspect, these simulations are not provided 
here. 
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Fig. 10: EV charging monetary costs at NE (converging point 
of Algorithm 1) against the number of EVs (/) when the 
individual payoff function (|5]) corresponds to: energy losses 
(a = 0 and = const.); transformer aging (a = 1, 

gf^ = const.); the charging monetary cost = 0 and 

French on/off peak fares (for 

and Scenario (s) {Ci = C = 16, Ui = a = 1, di = d = 30) 
are assumed. 



Fig. 11: This hgure correspond to the same setting as Fig. m 
except for the price model which is now the Epex Spot price 
model (i50ll. 

Price of Decentralization. Computing the PoD in general 
is a hard problem since it involves the maximization of 
the sum-payoff. Nonetheless, all the numerical results we 
obtained for various special cases allowed us to conhrm what 
Prop. IIV.4I suggests namely, that the PoD is close to zero. 
Concerning the case a = 1, it is seen from Pig. |6]that, under 
perfect forecasting, the transformer lifetime obtained thanks 
to Algorithm 1 is greater than 37.75 years while the no- 
EV upper bound is at 40 years. This shows that the PoD is 
necessarily less than — 1 ~ 6%. To elaborate further 

into this direction. Fig. [12] is provided in the case where 



Fig. 12: Price of Decentralization (PoD) in the simple sym¬ 
metric case L“° = 0, a = 0 (energy losses only), Oi = a = 1, 
di = d = T. The PoD is small (less than 20% for 7 = 4 EVs 
and C^=C = 2). 

a = 0. It represents the PoD against the number of EVs for 
a simplihed setting which allows exhaustive search for the 
centralized solution to be implemented and which is detailed 
on the hgure itself. Since quantities under use are discrete, 
there is a combinatorial effect which explains the different 
peaks. But the general tendency is that the PoD is relatively 
small and typically decreases with the number of EVs. All 
these observations fully support the relevance of distributed 
implementations of charging algorithms. 

VI. Concluding remarks 

In this section, we summarize a couple of key messages of 
this work and then provide possible extensions. 

One key feature of this work is to consider rectangular pro- 
hles for charging electric vehicles (EVs) and study in details 
the consequences of this choice. One important message is that 
rectangular prohles have the advantage to be more robust to 
errors on the non-EV demand forecast than existing charging 
schemes which rely on continuous power levels (e.g., valley- 
hlling solutions). As far as transformer lifetime is concerned, 
it is also seen that in the presence of forecasting noise, there 
exists an optimal power level at which the EVs should charge. 

It is seen that in the context of rectangular charging prohles, 
it is possible to capture in a simple manner the tradeoff be¬ 
tween transformer aging, energy losses, and charging monetary 
cost. As a rule of thumb for the engineer, it is observed that 
designing a good charging scheme in terms of aging will 
imply a good scheme in terms of energy losses whereas the 
converse does not hold in general (it may hold e.g., for a 
transformer with a very low thermal inertia). When electricity 
market prices and the non-EV demand are correlated e.g., 
by the importance of the evening peak, then a good scheme 
in terms of aging typically also performs well in terms of 
charging monetary cost. In more complex scenarios where 
this correlation is not available, the tradeoff may be analyzed 
through other simulations which are not provided here. 

























TRANSACTIONS ON SMART GRID 


12 


The distributed formulation of the charging problem in the 
considered context is shown to be fully relevant. Indeed, the 
proposed algorithm is shown to be scalable regarding conver¬ 
gence time, complexity, and required information. Concerning 
the specific issue of convergence, the ordinality potentiality 
property is proven to be available for key special cases 
and therefore guarantees convergence. Otherwise, in more 
general settings, simulations show that convergence is ensured 
with overwhelming probability. A possible downside of a 
distributed formulation of the problem is the potential loss 
of optimality. However, several strong arguments are provided 
to show that this issue is quite minor when the problem is 
formulated as in this work. 

A relevant extension is to use a dynamical approach which 
would take as a system state the transformer hot-spot temper¬ 
ature and charging policies as control policies. Since realistic 
models for the temperature evolution law are typically non¬ 
linear, the considered cost function is typically not quadratic, 
and charging constraints have to be considered, the underlying 
control problem seems to be non-trivial. Using a dynamical 
approach is also relevant to account for the memory effect 
in electricity prices. Additionally, the problem of robustness 
against forecasting errors has to be accounted for. Therefore a 
quite complete mathematical model would be to account both 
for memory effects (given by the transfomer and the prices) 
and uncertainty, which would lead to using the advanced tools 
of stochastic control 15^ , 15^ and games. 


Appendix A 

Proof of the potential property of the charging 
GAME (Prop. IV. 3) 

Proof: We start with the case of assumption (a) i.e., Vi S 
I, Wi{si) = W. In this case, we have the following sequence 
of statements; 

1 ) the game in which the payoff function is the same for all 
EVs and given by M^(s) = m^(s) = aAt [L* (s)) + 

tew 

(1 — a)J {Lt (s)) is an (exact) potential game (it 
is a "team game" more precisely); 

2) the game in which the payoff function of EV i is 
uj{s) = gf^ (si) is an (exact) potential game because 
the payoff function of each player only depends on its 
own strategy; 

3) the game in which the payoff function of EV i is 
u^{s) = ul{s) + uf{s) is an (exact) potential game 
because the payoff functions of the players are the sum 
of two payoff functions for which the game is an (exact) 
potential game; 

4) the game in which the payoff function of EV i is 
uj{s) = —/i{itf(s)) is an ordinal potential game 
because composing the payoff functions of a potential 
game by strictly decreasing functions leads to an ordinal 
potential gamqj. 


We now treat case (b). We show that 

nt(s) 

^(6) (s) = - X! + (1 - a)J {Lt{vt)) 

teT vt=o 

(15) 

iei 


is an ordinal potential for the charging game. Suppose that 
player i deviates from Si to s^. Assume, for example, that 
this deviation is rational and brings him a payoff increase of 
-/. ( 5 ,™ (s', s_.) + gr (s')) + h {gf^ (s) + gf (s.)) > 0 
We have the following relations 


fi (s) + gr (si)) 


- /* {g?^ {s'i, s-i) + g^ {s'i)) > 0 
)+5r(so-<?r(s')>o (16) 


t=Si t=s[ 

+5r(s.)-5r(s')>o (17) 

O (s') _ (s) > 0 , (18) 

where (L* (s)) = aAt {Lt (s)) + (1 - a)J {Lt (s)) to 
simplify the notations. 

The first equivalence (fThl) comes from the strict mono¬ 
tonicity of fi, while (fTTl) and (fTSl l follow by observing 
that the only terms changing in because of player Es 
deviation are nsi,nsi+i, ...Us^+Ci-i (decreasing by one), and 
nst,nst+i, ■■■,ns'+Ci-i (increasing by one), which modifies 
the corresponding sums in (fTsT i by, respectively, subtracting a 
term for t = Si, Si + 1, ...Si + Ci — 1, and adding a term for 
t = s', s'-fl, ...s'-l-Ci —1. It therefore gives the equivalence in 
the Ordinal Potential Game definition (Def. |2]i and concludes 
the proof. ■ 


Appendix B 

Proof of the PoD in a special case 

Proof: We consider scenario (b) for the proof. The proof 
for a = 0 is obtained directly from this one observing that 
setting q; = 0 leads to a DN cost without inertia. Again, let 
g^^ {Lt (s)) = aAt {Lt (s)) -f (1 - a)J {Lt (s)) to simplify 
the notations. By following the steps in the proof of App. 
it can be shown that the nonatomic version of the charging 
game is also an OPG, with a potential function candidate 


root 

= g^^ {Lt{vt)) dvt (19) 

t^'fJvt=0 

where Lt{vt) = Lf'° + pvt ■ Note that there is no individual 
term related to gf^ in ( fT9] l given the assumption made in Prop. 
121 

Observe now that without the constraints xi < ■ ■ ■ < xc 
and xt-c +1 > •' • > xt, which ensure that no EV 
starts charging before f = 1 or finishes after t = T, the 
minimization problems of both —w on the set 

X = {a; £ Xt = C} are standard valley¬ 

filling problem^ with a strictly convex cost function, here 


^mina, +pxt) subject to a: > 0, xt = C, with / the cost 

^Furthermore, when fi = Id, the charging game is an exact potential game. function. 




TRANSACTIONS ON SMART GRID 


13 


X I—5^^ {Lt{v)) di; in the case of is strictly 

increasing) and x i— {Lt{x)) in the case of w. It is known 
that the solution of this problem is unique and independent of 
the strictly convex cost function (see e.g., Thm.2 in ITSl ). Let 
X* denote this solution, x* has the following "valley-filling" 
structure; 

yt, xt>o^Lr + pi£; = L* 
yt, xt = 0^ Lf° >L* ’ ^ 

where L* denotes the "valley level" of the solution. This 
implies here that the solutions of the minimization problems 
of both and re on X coincide. 

It remains to show that under the assumptions Lf'° > ■ ■ ■ > 
Uq° and < • ■ • < L™, x* is still the solution of 

both problems on X (instead of X), which is equivalent to 
X* G X. In fact, it is easy to see that x^ < ■ ■ ■ < x^ and 
x’^-C+i — ■ ■ ■ ^ necessarily hold. Suppose indeed that 
the converse is true e.g., x^ > This implies L™ + pxl > 

Lf°+px 2 which in turn implies = 0 because of ( l20l i, which 
is contradictory. Then, x* is the (unique) solution of both 
optimization problems on X. Thus, the (unique) NE profile is 
maximizing w on X and PoD = 0. ■ 
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